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We report an intrinsic absorber with decreasing outflow velocity in the Seyfert 
1 galaxy NGC 3783. This is the first detection of a change in radial velocity in 
an outflow associated with a Seyfert galaxy. These results are based on mea- 
surements from 18 observations with the Space Telescope Imaging Spectrograph 
aboard the Hubble Space Telescope, obtained between 2000 February and 2002 
January. In two intervals separated by ~ 13 and 9 months, the absorption lines 
in the kinematic component with highest outflow velocity exhibited mean red- 
ward velocity shifts of ~ 35 and 55 km s -1 , respectively. The rate of velocity 
decrease was 2.2 ± 0.6 times more rapid in the second interval. No variations in 
absorption velocities were detected in the other kinematic components. We ex- 
plore potential interpretations of the observed velocity shifts: radial deceleration 
of the UV absorber due to a change in either the speed or direction of motion of 
the outflow, and the evolution of a continuous flow across our line of sight to the 
emission source. 

Subject headings: galaxies: individual (NGC 3783) — galaxies: active — galaxies: 
Seyfert — ultraviolet: galaxies 

1. Introduction 

Recent surveys with the Hubble Space Telescope (HST) and the Far Ultraviolet Spectro- 
scopic Explorer (FUSE) have shown intrinsic UV absorption is common in Seyfert 1 galaxies, 
appearing in over half of the objects observed to date (Crenshaw et al. 1999; Kriss 2002). 
In most cases, the absorption is blueshifted with respect to the systemic velocity of the host 
galaxy, indicating radial outflow from the active galactic nucleus (AGN). There is a one- 
to-one correspondence between the detection of intrinsic UV absorption and X-ray "warm 
absorption" in Seyfert galaxies, indicating they are closely connected (Crenshaw et al. 1999). 
The mass outflow revealed by the blueshifted absorption may be a major component of the 
energetics and circumnuclear environment in AGN (e.g., Mathur et al. 1995; Reynolds 1997). 
Currently, the mechanism driving this mass outflow in AGN is not well understood. 

Variability of the intrinsic UV absorption appears to be common in Seyferts. In the few 
objects that have multiple UV spectra at sufficiently high spectral resolution for the accurate 
comparison of column densities, significant variations are detected in all cases (Maran et 
al. 1996; Crenshaw et al. 2000a; Kraemer et al. 2002; Crenshaw et al. 2003). In some 
extreme examples, entire kinematic components have appeared and disappeared on roughly 
yearly timescales. Strong evidence for column density variations ascribable both to ionization 
changes in response to the AGN continuum (Kraemer et al. 2001a, 2002) and to changes in 
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total absorption columns in our line of sight to the AGN, i.e., due to transverse motion 
(Kraemer et al. 2001b; Crenshaw et al. 2003), have been found. 

One somewhat surprising finding has been the lack of variability in the absorption 
kinematics. Although multi-epoch data are sparse, no detections of radial velocity variations 
have been reported for any intrinsic absorbers in Seyfert galaxies, revealing stability on 
timescales of at least several years in some cases (e.g., Weymann et al. 1997; Kraemer et al. 
2001b). Many of the extreme counterparts to Seyfert absorbers, the broad absorption line 
(BAL) systems seen in more luminous AGN, have also shown this long-term stability (e.g., 
Barlow 1993; Weymann 2002). Limits on radial acceleration in BALs have been used to 
assess models of the dynamical forces associated with their mass outflow (Foltz et al. 1987; 
Weymann et al. 1997). Rare exceptions to this stability have been reported by Hamann 
et al. (1997) and Vilkoviskij & Irwin (2001), in which centroid radial velocity shifts were 
detected in relatively narrow absorption line systems in QSOs. 

The bright Seyfert 1 galaxy NGC 3783 has exhibited dramatic variability in its UV 
absorption spectrum. In only four observations with the Goddard High Resolution Spec- 
trograph (GHRS) and the Space Telescope Imaging Spectrograph (STIS) obtained approxi- 
mately 1-5 years apart, new kinematic components appeared in three separate cases (Maran 
et al. 1996; Kraemer et al. 2001b). No correlation was found between the column densities 
and continuum flux, indicating these large-scale changes were primarily due to transverse 
motion of the absorbers (Kraemer et al. 2001b). In part because of this extreme variability, 
NGC 3783 was selected for an intensive multiwavelength monitoring campaign with HST, 
FUSE, and Chandra. In this study, we concentrate on variability of the absorption kine- 
matics based on analysis of the STIS spectra obtained during our monitoring campaign. 
Analyses of the averaged X-ray and UV (STIS and FUSE) spectra are presented in Kaspi et 
al. (2002, hereafter Paper I) and Gabel et al. (2003, hereafter Paper II), respectively. Full 
treatment of the variability of the UV and X-ray absorption will be given in subsequent 
papers. 



2. Observations and Results 

A total of 18 HST '/STIS observations of the nucleus of NGC 3783 were obtained between 
2000 February 27 and 2002 January 6 using the E140M medium resolution echelle grating. 
This grating spans 1150 - 1730 A with a spectral resolving power of R ~ 40,000, sampling 
the Lya, N V AA1238.8,1242.8, C IV AA1548.2, 1550.8, and Si IV AA1393.8,1402.8 resonance 
lines at a resolution of ~8 km s -1 . Each visit consisted of two orbits and ~4.9 ks of exposure 
time. A log of the observations and full description of the data reduction are given in 
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Paper II. 

A detailed discussion of the UV absorption spectrum is presented in Kraemer et al. 
(2001b) and Paper II. Here we review some features that are important for the subsequent 
analysis. During our monitoring campaign, strong absorption was present in Lja, N V, and 
C IV in three kinematic components, at centroid radial velocities of v r —1350, —725, and 
—550 km s -1 relative to the systemic redshift of NGC 3783 (z = 0.00976; de Vaucouleurs et 
al. 1991). For historical reasons, we refer to these as components 1, 3, and 2, respectively. 
Weak absorption was tentatively detected in a fourth component at —1050 km s _1 . The Si IV 
doublet only appears in component 1. Metastable C III was detected in component 1 in the 
averaged spectrum, indicating a high density in this absorber (n e ~ 10 9 cm -3 ; Paper II). 
Absorption in these components first appeared in the following epochs: component 2 in a 
1994 January 16 GHRS spectrum of C IV, component 1 in a 1995 April 11 GHRS spectrum 
of C IV, and components 3 and 4 in the first STIS spectrum (2000 February 27). 

We fit the intrinsic (i.e., unabsorbed) continuum plus emission-line fluxes in each STIS 
spectrum following the procedure described in Paper II. Each spectrum was normalized by 
dividing by its fit. We compared the interstellar absorption lines in each spectrum to ensure 
that there were no wavelength offsets. Specifically, we measured the centroid wavelengths 
of the Galactic Si IV A1393.8, C IV A1548.2, and Fe II A1608.5 lines in all spectra. For 
each interstellar line, the standard deviation of the centroids (~ 2 - 4 km s _1 ) is less than a 
resolution element, indicating accurate relative wavelength solutions for these observations. 
Centroid radial velocities were then measured for each intrinsic absorption line. 

In Figure la, the N V, C IV, and Si IV component 1 centroid velocities measured in 
all STIS spectra are plotted, clearly showing the radial velocity has decreased during the 
monitoring campaign. We identify three epochs exhibiting marked velocity variations: the 
first observation obtained 2000 February 27, the intensive phase of the campaign consisting 
of 13 observations obtained between 2001 February 25 - April 23, and the final observation 
obtained 2002 January 6. Figure lb show only the intensive monitoring observations (epoch 
2). Careful inspection of Figure lb reveals no correlations in the centroids measured for the 
three ions; the fluctuations are likely due to the measurement uncertainties (see discussion 
below) . 

Table 1 lists the component 1 velocity measurements for each of the three epochs iden- 
tified above. Centroids and uncertainties quoted for epoch 2 are the means and standard 
deviations of the 13 observations obtained during the intensive phase of the campaign. These 
observations also represent a good statistical sample for determining the random uncertain- 
ties associated with measurements in the individual spectra. Thus, we adopted the standard 
deviations for the 1 a errors in epochs 1 and 3 as well, with the caveats that (1) if there are 
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real variations in the centroid velocity during the intensive phase of the campaign these are 
upper limits on the random uncertainties and (2) this does not account for any systematic 
errors that may be associated with the measurements. The radial decelerations listed in 
Table 1, a r , assume a constant decrease in velocity between epochs. Uncertainties in these 
values were derived by propagating the errors associated with the radial velocities. The 
average deceleration between epochs 2 and 3 was 2.2 ± 0.6 times greater than in the interval 
between epochs 1 and 2. N V shows some evidence for smaller velocity shifts than the other 
lines; however, formally, the decelerations are equal within the 1 a uncertainties. No radial 
velocity variations were detected in kinematic components 2-4 within the uncertainties 
in the measurements (Av r ~ 10-20 km s^ 1 ). Interestingly, only the component with the 
highest outflow velocity is observed to be decelerating. 

While a full treatment of the other absorption properties, i.e., absorption strengths, 
widths, and covering factors, will be given in a later study (Gabel et al., in preparation), we 
briefly describe results here that have implications for the interpretation of the radial velocity 
shift. First, we note that although there were some variations in the absorption strengths 
and widths during the monitoring campaign, there were no large-scale changes in the three 
epochs that would indicate we are seeing different components appearing and disappearing 
to mimic a shift in velocity (although we may be seeing different regions of a continuous 
flow moving across our line of sight, as discussed in §3.3). This is evident from the overall 
similarity in the absorption spectra in Figure 2, which shows the profiles for each line in 
each of the three epochs. There were moderate variations in the component 1 line widths 
during the STIS observations, but these variations were complex; each line varied differently 
and there were no apparent correlations with radial velocity (FWHM are listed in Table 1). 
For example, Figure 2 shows the mean C IV width in epoch 2 (FWHM =150 km s _1 ) is 
narrower than in epochs 1 and 3 (FWHM =203 and 194 ±15 km s _1 , respectively), however 
the width varied between 115 - 190 km s _1 during the epoch 2 observations, on much shorter 
timescales than the radial velocity shifts. Furthermore, Table 1 shows the N V widths were 
not correlated with C IV. On the other hand, the C IV profile in the 1995 GHRS spectrum 
is much broader (FHWM =500 ±50 km s _1 ) than during the STIS observations, as seen in 
Figure 1 in Kraemer et al. (2001b). Given the five year separation between the GHRS and 
STIS observations and the dramatic difference in the profiles, we cannot rule out that the 
absorption in the GHRS spectrum is unrelated to that observed with STIS; however, if it is, it 
indicates a significant evolution of the absorber's structure over this interval. There were no 
detectable changes in the line-of-sight covering factors in the component 1 absorption lines. 
This provides constraints on the evolution of the absorption/emission geometry between the 
three epochs. Finally, we note there is evidence that the UV absorption in component 1 
is comprised of at least two subcomponents with different ionization structures. This was 
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initially suggested by Kraemer et al. (2001b) from photoionization modeling of the first STIS 
spectrum. Also, the Si IV covering factor in component 1 was found to be smaller than that 
of at least H I and N V (Paper II). Since all lines exhibit the velocity shift, this indicates 
both subcomponents are decreasing in outflow velocity. 

3. Interpretation of Decreasing Radial Velocity 

3.1. Deceleration of the Outflow? 

One possible explanation for the observed decrease in radial velocity is that the absorber 
is undergoing a bulk radial deceleration. We first consider gravity as the source of the 
required inward radial force. The simplest case, gravitational deceleration of a ballistically 
ejected absorber, is excluded by the increased deceleration detected in the second interval 
(see Table 1), since the radial distance between the absorber and central mass has increased 
between epochs (by ~3xl0 15 cm). This requires the presence of an outward radial force 
that has decreased relative to gravity between the two intervals. Since the other kinematic 
components in NGC 3783 are not decelerating, this adds the further requirement that the 
nature of the radial forces for these components differs from those of component 1, e.g., 
they may have reached sufficient distances from the central mass that they are near terminal 
velocity. Using the average deceleration over the second interval listed in Table 1 and the 
black hole mass for NGC 3783 (M BH =9xl0 6 M ; Onken & Peterson 2002), the absorber 
would be < 9xl0 16 cm from the central mass if gravity is the decelerating force (the upper 
limit is for a negligible outward force during this interval). 

An alternate decelerating mechanism could be an interaction of the absorber with other 
gas in the vicinity of the AGN. This has been invoked to explain observations of deceleration 
in the narrow line region (NLR) gas in several Seyfert galaxies because the onset of decel- 
eration in these objects is at distances well beyond the gravitational influence of the central 
mass (Crenshaw & Kraemer 2000; Crenshaw et al. 2000a; Ruiz et al. 2001). In the context 
of intrinsic absorbers, this could be the signature of an interaction between different phases 
in the outflow, as predicted in many dynamical models (e.g., Begelman et al. 1991; Krolik 
& Kriss 2001; Everett 2003). The hydrodynamic drag force is proportional to Av 2 rii/N a , 
where Av is the difference in velocity between the absorber and the ambient medium, rij is 
the density of the ambient medium, and N a is the column density of the absorber (Begelman 
et al. 1991). In this framework, a twofold increase in rii/N a is required in the second interval 
to account for the observed increased deceleration. However, a potential problem for this 
model is the absence of any signatures of the interaction in the absorption spectrum. The 
deceleration of the NGC 3783 absorber is much greater than that observed in the NLR gas, 
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requiring a drag force (per unit mass) that is four orders of magnitude greater (F^/m ~ 0.05 
cm s~ 2 ). Changes in the ionization structure appear in the decelerating region of the NLR 
gas, presumably due to shocks (Kraemer & Crenshaw 2000; Cecil et al. 2002), whereas no 
such changes are detected in the absorption spectrum in NGC 3783, as discussed above. 

A third possible mechanism is deceleration of the outflow due to mass loading. In 
this scenario, additional mass is swept up by a momentum conserving flow, resulting in a 
decreased outflow velocity. This is believed to be the dominant effect in the deceleration 
of supernova remnants (e.g., Dyson et al. 2002, and references therein). As an example of 
how an outflow's mass could be increased, we note the dynamical models of Everett (2003) 
predict material from the accretion disk is continually lifted by a disk-driven wind and added 
to the flow. Since the deceleration will depend on the rate of increase in the flow's mass, the 
increased deceleration observed in the second interval requires that mass was accumulated 
more rapidly than during the first interval. 



3.2. Directional Shift in the Outflow? 

Another interpretation of the observed decrease in blueshifted radial velocity is that, 
rather than a bulk radial deceleration as described above, we are seeing a change in the 
radial component of the velocity vector with respect to our line of sight to the AGN. Figure 3 
illustrates this schematically, depicting the absorber as a cloud moving along a curved path 
across our sight line to the background emission source. In this simple model, the velocity of 
the cloud, represented with thick filled arrows in Figure 3, has a constant amplitude along 
the trajectory, but the component that is directed radially in our line of sight (open arrows 
in the x r direction) decreases due to the changing direction of the motion. The absorber 
occults ~30% of the BLR (Paper II), which is demarcated with dashed lines in Figure 3. 

The transverse component of motion of the absorber (i.e., along xt in Figure 3) is 
an integral part of this model, leading to the motion of the absorber into and eventually 
out of our view. A rough lower limit on the transverse velocity vt can be derived using the 
transverse size of the absorber (lxlO 16 cm; Paper II) and an upper limit on the time required 
for the absorber to move into our line of sight. As a conservative estimate of the latter, we 
use the interval between the 1994 January GHRS observation and the first STIS observation, 
since it is uncertain whether the absorption in the 1995 April GHRS spectrum is from the 
same absorber as that in the STIS data (see discussion in §2). The result is vt > 542 km s _1 , 
which, in combination with the transverse size of the BLR (1.9 xlO 16 cm; Onken & Peterson 
2002), gives a rough upper limit on the total time that some portion of the absorber will 
occult the BLR, t < 17 yrs. This is strictly an upper limit since (a) we have only an 
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upper limit on the interval over which the absorption appeared and (b) Vt is increasing, as 
illustrated in Figure 3. Since the continuum source is certainly much smaller than the BLR, 
at some point in the absorber's path it will no longer occult the continuum source but will 
still occult the same fraction of the BLR (see Paper II for a full discussion of the covering 
factors). Thus, this model predicts the covering factors, and therefore absorption depths, 
could change substantially at a later epoch, depending on the specific geometry of the BLR 
and continuum sources. Since the covering factor was not observed to change between the 
three epochs in our STIS campaign as described in §2, this implies the absorber has not yet 
moved out of our line of sight to the continuum source. We note that the rate of change in 
radial velocity will be determined by the curvature of the absorber's trajectory and will not 
necessarily be constant. This may provide a natural explanation for the more rapid decrease 
in radial velocity observed in the second interval of the observations. 

A number of mechanisms could produce the curved motion depicted in Figure 3. For 
example, models positing the absorber is driven off the accretion disk by radiation pressure 
from the central source (e.g., Murray et al. 1995; Proga et al. 2000) or by centrifugal forces 
along magnetic field lines in the disk (e.g., Blandford & Payne 1982; Bottorff et al. 2000; 
Everett 2003) all predict the outflow will follow a curved trajectory above the disk, exhibiting 
a changing angle 9 with respect to the perpendicular to the disk. Additionally, these models 
predict the outflow's path will be helical due to rotation about the central mass, and thus 
its motion will be curved in the plane parallel to the disk. 



3.3. Evolution of a Continuous Flow? 

If the component 1 absorber is a continuous outflow, the centroid velocity could vary due 
to the motion of different regions of the flow across our line of sight to the emission source. In 
this scenario, the absorber could be either a continuous stream of gas or a continuous "spray" 
of small clouds. An example is shown in Figure 4, in which a collimated flow, dileneated with 
thick solid lines, intersects our line of sight to the emission source obliquely. The mean angle 
of intersection in the flow, 9, determines the centroid radial velocity. Due to its finite opening 
angle, A9, the flow spans a range of projected radial velocities, giving the observed width 
of the absorption feature, v x (cos(9 + Ad/ 2) — cos{9 — A9/2)). Two epochs are depicted 
in Figure 4, showing how changes in the flow geometry lead to variations in the absorption 
kinematics. Centroid velocity variations could be produced by either an overall change in 
the flow's direction or inhomogeneities in the flow within a given opening angle that result 
in a change in 9. In constrast to the models discussed in the previous sections, the observed 
velocity shift is due to the motion of different gas elements, with different projected radial 
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velocities, into our line of sight rather than the change in speed or direction of motion of any 
given gas element. 

4. Conclusions 

STIS spectra obtained as part of our monitoring campaign of NGC 3783 reveal a de- 
creasing blueshifted radial velocity in one of its intrinsic absorption components. This is the 
first detection of a variable radial velocity in an outflow associated with a Seyfert galaxy. 
In two intervals separated by ~ 13 and 9 months, the absorption lines in the kinematic 
component with highest outflow velocity exhibited mean redward velocity shifts of ~ 35 
and 55 km s _1 , respectively. The rate of decrease in the second interval was twice as large 
as in the first. The absorption widths in this component varied moderately during these 
observations. These variations were complex, fluctuating on short timescales compared to 
the velocity shifts and showing no apparent correlation with the variations in radial velocity. 
No variations in absorption velocities were detected in the other kinematic components. 

We have explored various interpretations of the observed velocity shift. The similarity 
in spectra at all epochs indicates it is not due to changes in the relative strengths of multiple 
subcomponents, and it thus represents a bona fide decrease in radial velocity. If it represents 
bulk deceleration due to a decrease in radial speed, we consider gravity, interaction with an 
ambient medium, and mass loading as possible deceleration mechanisms and derive special 
requirements on the nature of the radial forces. Alternatively, it could be due to a directional 
shift in the motion of the absorber with respect to our line of sight to the background emission 
sources, or a result of different regions in a continuous flow crossing our sightline. 

J. R. G., D. M. C, and S. B. K. acknowledge support from NASA grant HST-GO- 
08606. 13- A. W.N.B acknowledges CXC grant GO 1-2 103 and NASA LTSA grant NAG5- 
13035. F.H. acknowledges NSF grant AST 99-84040. We thank the anonymous referee for 
very helpful comments. 
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Table 1. Variable Kinematics in Component 1 



Epoch 


Julian Date 




N V 






C IV 






Si IV a 










FWHM 






FWHM 








(2,450,000 +) 


(km s -1 ) 


(10~ 6 km s" 2 ) 


(km s -1 ) 


(km s~ 1 ) 


(10~ 6 km s~ 2 ) 


(km s _1 ) 


(km s -1 ) 


(10~ 6 km s" 


1 


1602 


-1334±4 




212±15 


-1352±8 




203±15 


-1374±13 




2 d 


1966-2023 


-1307±4 


-0.8±0.2 


140-212 


-1318±8 


-1.0±0.3 


115-190 


-1333±13 


-1.2±0.5 


3 


2281 


-1261±4 


-1.8±0.3 


150±15 


-1256±8 


-2.5±0.5 


194±15 


-1280±13 


-2.2±0.8 



a Average values for the two lines of the Si IV doublet. 

b Ccntroid velocity of component 1 relative to the systemic velocity. 

c Average radial deceleration since previous epoch. 

d v r and a r are mean values and FWHM are range of values measured during intensive phase of monitoring. 
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Fig. 1. — Component 1 radial velocities measured in the STIS spectra. Centroid velocities 
and error bars are plotted for N V A1238.8 (top), C IV A1548.2 (middle) and the average of 
the two Si IV doublet lines (bottom) as a function of Julian date. Measurements from all 
observations are shown in (la); only observations obtained during the intensive monitoring 
are shown in (lb). 
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Fig. 2. — Normalized N V, C IV, and Si IV absorption profiles for three epochs showing 
the decrease in the component 1 radial velocity. Profiles are plotted in velocity space with 
respect to the short wavelength lines of the doublets in the rest frame of NGC 3783. Epoch 
1 is the first observation (dashed line), epoch 2 is the average of 13 observations obtained in 
the intensive phase of the campaign (solid line), and epoch 3 is the final observation (dotted 
line). All intrinsic absorption features are marked with tick marks above the spectra and 
numbered following the conventions described in the text: lower, thick marks identify the 
short wavelength doublet members; upper, thin marks identify the long wavelength members. 
The Si IV profile is the mean of the two doublet lines. Interstellar absorption features are 
marked with an "x" below each spectrum. Epochs 1 and 3 have been smoothed for clarity, 
leading to shallower interstellar lines compared to the composite spectrum of epoch 2. 
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Fig. 3. — Schematic representation of a simple model in which the curved path of an ab- 
sorption cloud leads to the observed decrease in radial velocity. The dashed lines outline 
the transverse size of the BLR projected on the sky. The coordinate system, relative to the 
observer, is shown at the top right. The absorber is depicted as a shaded circle that occults 
~30% of the BLR and moves along a curved path at constant velocity (thick, filled arrows). 
The radial component of the velocity (thin, open arrows) is seen to decrease at different 
points on the clouds trajectory. 




Fig. 4. — Schematic demonstrating how kinematic variability could arise in a continuous 
flow. Two epochs are shown, ti and t2- In each epoch, a collimated outflow, delineated 
with thick solid lines, intersects the observers line of sight (dashed lines) at an oblique angle. 
Due to the opening angle of the flow, a range of angles, A9, is spanned in the line of sight 
giving the absorption line width. The mean intersection angle, 9, determines the centroid. 
The resulting variability in absorption profiles due to the geometrical differences in different 
regions of the flow is shown at the top. 



